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The relativistic flux tube model is an effective description of confined quarks and gluons in which 
the confining interaction is carried by the flux tube, a Nambu-Goto string. We first show that the 
relativistic flux tube model can be applied to glueballs seen as bound states of transverse constituent 
gluons. After a comparison of that approach with usual spinless Salpeter Hamiltonians, we compute 
glueball masses and decay widths. Comments about the 7;-r;'-pseudosclar glueball problem, the 
glueball-Pomeron conjecture, and finite-temperature effects are finally given. We also point out the 
existence of a duality between open- and closed-flux tube models of glueballs. 

PACS numbers: 12.39.Mk 



I. INTRODUCTION 



Quantum chromodynamics (QCD) allows the existence of pure gauge states called glueballs. Although much 
effort has been devoted to find a clear experimental glueball signal, no unambiguous candidate is known so far 
Nevertheless, glueballs are intensively studied within various theoretical frameworks: lattice QCD [1, 0, 0], Coulomb 
gauge QCD ;5|, AdS/QCD 6], QCD sum rules 0, potential models [1, i, Ell, d , etc. We refer the interested 
reader to Ref. |13| for a review on theoretical glueball physics. Nowadays, the glueball mass spectra computed in 
lattice QCD have a special status since they are often taken as input to fit the parameters of the existing models. 
Indeed, clear experimental data related to glueballs are still missing. 

The present work is devoted to an effective approach of QCD that has poorly been applied to glueballs so far: the 
relativistic flux tube model [13, IH, [l^- Actually, in the meson case, it follows from nonperturbative QCD in the 
quenched approximation that the confining interaction between a quark and an antiquark pair is due at the dominant 
order to a Nambu-Goto string linking them and carrying both energy and angular momentum p^ . This agrees 
with lattice QCD simulations showing that the chromoelectric field between a static quark-antiquark pair is roughly 
constant on a straight line joining these two particles (see for example the early study [T3|, and the more recent works 
[isl. fl9l|). Calculations within the dual Ginzburg-Landau theory also support the existence of a flux tube [lO]- The 
relativistic flux tube model is not a potential approach since the interaction is mediated by a dynamical object. In 
particular, not only light and heavy mesons can be described in that framework [2ll [2^ . but also hybrid mesons 
when the flux tube is in an excited state [11]. Here we extend the relativistic flux tube model to glueballs, assuming 
that these exotic hadrons can be described as bound states of transverse constituent gluons. We point out that using 
longitudinal constituent gluons leads to spurious states with respect to those observed in lattice QCD, and to a mass 
spectrum that cannot agree with the lattice data in the C = — sector [13, [l^ ■ 

After having recalled how to build quantum states of transverse gluons in section |TT1 we present the relativistic flux 
tube model and its generalization to gluonic bound states in section Hill We also study the Regge trajectories, that 
are a consequence of the relativistic flux tube model in the ultrarelativistic limit, in section Hvl Then we compute 
the glueball mass spectrum in section |V] and estimate the glueball decay widths in section IVII Finally, we discuss 
about the gluonic content of the rj-rj' resonances, the glueball-Pomeron conjecture, and finite-temperature effects in 
sections IVID I Villi and IIXI Our results are summarized in section |Xl 
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II. BOUND STATES OF TRANSVERSE GLUONS 



A. Generalities 

The main goal of this work is to build an effective model of glueball based on the relativistic flux model in which 
the constituent gluons are transverse, that is with helicity-1. As we argued in (T^. [23|. the helicity formalism of Jacob 
and Wick [2^ can be successfully applied to describe a bound state of two such gluons. It is thus important to show 
here how to build glueball helicity states before using them as basis states for computations within the relativistic 
flux tube model. 

Let us recall the main points of the helicity formalism, introduced in [25} . We first introduce A)) = a\^{p ) |0) 
the quantum state of a particle with momentum p, spin s, and helicity A. If the particle is transverse (with a helicity 
degree of freedom), only A = ±s is allowed, while the (2s + 1) projections from —s to +s are allowed if the particle 
has a spin degree of freedom. Then it can be deduced from [25| that the quantum state 
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represents a two-particle state in the rest frame of the system which is also an eigenstatc of the total spin J, i.e. 

— J {J + 1) and Jz = M. In the above definition, R{a, (3, 7) is the rotation operator of Euler angles {a, /3, 7} and 
V'lf ^{a, (3, 7) are the Wigner Z?-matrices. The coordinates {6, 0} are the polar angles of p. It is worth recalling the 
inequality 

J>|Ai-A2|, (2) 
coming from usual spin-coupling rules. Sill by applying the formalism of [1^, it can be checked that the state 

|Ai,A2; J^,M,e) = -i= [|Ai, A2; J, M) -I- e j-Ai, -A2; J, M)] (3) 

is also an eigenstate of the parity operator with the eigenvalue 

P = e^,Tj2{-iy-'^-'\ (4) 

rji and Si denoting the intrinsic parity and spin of particle i respectively. Let us note that only the relative sign of Ai 
and A2 is relevant. 

Either particles with spin or helicity can be described within this formalism. When both particles have a spin degree 
of freedom, the helicity basis, spanned by the helicity states (O, is equivalent to a standard p'^+^Lj) basis up to an 
orthogonal transformation [2y|. When at least one of the particles is transverse, both basis are no longer equivalent 
but the helicity states can still be expressed as particular linear combinations of p'^+^Lj) states [H, [25|. It is worth 
mentioning that only the spin-angular part of the helicity states is fixed by the above geometrical construction. The 
radial part of the system needs a particular dynamical model to be determined. 

Although we will mainly focus on two-gluon states in the present work, it is worth mentioning that three-body 
helicity states can also be built following a two-by-two recoupling scheme [1^ [l^. The idea is to write first the 
two-body helicity state ([3]) corresponding to, say, particles 2 and 3 in their rest frame. Then, the Lorentz-boosted 
helicity state of the (2,3) cluster can be coupled to particle 1 with the same technique. The boost is needed for the 
three-body helicity state to be expressed in the three-body rest frame, not in the rest frame of the cluster. That 
recoupling procedure can in principle be generalized to A/'-body systems. 

B. Two-gluon states 

Transverse gluons are such that A^ = ±1. Since e = ±1, there are thus four independent helicity states of the 
form ([3]), and the Pauli principle must be satisfied. The color wave function of a two-gluon system, [8, S]-*", is symmetric 
and enforces a positive charge conjugation; one has thus to ask the two-gluon states to be totally symmetric. The 
symmetrization leads in general to selection rules on the total spin of the system [2^. We have shown in Ref. 
that, in the case of two gluons, the symmetrized helicity states read 

|^+;(2A;)+) = 1 1, 1; (2A:)+, Af, l) , \D+; {2k + 2)+) = 1 1, -1; (2fc + 2) + , Af, l) , (5) 

\S-;(2k)-) = |l,l;(2/c)-,A'/,-l), |D_;(2fc + 3)+) = 1 1, -1; (2fc + 3) + , Af , -l) , 



with k E N. The S (£')-states will be referred to as helicity singlets (doublets) in the following. It appears that there 
are no J^^ = two-gluon states, in agreement with Yang's theorem forbidding the decay of vector mesons in two 
photons [11]. This fact is also consistent with the absence of low- lying 1^+ glueballs in the lattice QCD results 

Once the glueball helicity states ([S]) are known, it is possible to compute matrix elements. In particular, we have 
checked in Ref. that none of the usual operators (identity, spin-orbit, tensor, . . . ) induce couplings between the 
S+ and £>+ states, that possess the same J^'-^ quantum numbers. Moreover, the average value of the square orbital 
angular momentum takes the simple form 

<Ji2^== J(J+l) + 2AiA2. (6) 

It can be expected that (^-^^^ roughly sets the energy scale of a given state: The more rotational energy is contained 

in a glueball, the heavier the state should be. This suggests the following mass or dering Q: 0±+, 2++, 2"+, 3++, 
etc., in agreement with the lightest glueball states observed in lattice QCD [2]. Two is indeed the minimal number of 
gluons needed to make a color singlet, and it is rather intuitive to see the lightest glueballs as two-gluon bound states. 



C. Many-gluon states 

The situation is more complicated for thrce-gluon glueballs. First of all, two color singlets are possible: One 
symmetric (leading to C = — ) and one antisymmetric (leading to C = +). The corresponding color wave functions are 
[[8, 8] ° , 8]-'-= and [[8, 8] " , 8]-'-'' respectively. The low-lying C = — glueballs are well-known from lattice computations. 
We gave in Ref. [2J] several arguments favoring a three-gluon interpretation of these states, the most obvious ones 
being that negative charge conjugations cannot be reached by a two-gluon bound state, and that the lightest C = — 
glueballs are heavier than the C = + ones. Moreover, the lowest-lying bound states of three gluons should be J*^ — 1^ 
and 3~ ones, as observed in lattice QCD [24j. Nevertheless, writing completely symmetrized three-gluon states in the 
helicity formalism is a complex problem that would deserve a separate paper; it is beyond the scope of the present 
work. 

It is worth saying that building a (pseudo)scalar state from a totally symmetric three-gluon system might be 
problematic. This case is indeed formally identical to a three-photon system. It has been previously shown that a 
planar three-photon system where the photons have equal energies and momenta separated by 120° from one another 
cannot be in a (pseudo)scalar state [23|- Such a situation is analogous to low- lying three-gluon states; no light 
glueball is indeed observed. However, asymmetric three-photon configurations, that could correspond to highly 
excited glueballs, can be in such a state [3(|. Consequently, if there is a three-gluon glueball with J'^ = 0~, it must 
be particularly heavy. In agreement with that picture, the heaviest-known glueball on the lattice is a 0^ state. But 
we argued in Ref. [24j that, rather than a highly excited three-gluon state, the 0^ glueball could also be a low- lying 

four-gluon state with the nontrivial color wave function [[8, 8] " , [8, 8] "]^ imposing the mixed symmetry 



to 



the spin-space wave function. Notice that C — + glueballs in the mass range of the C = — ones are currently not 
known from lattice computations since one- and two-glueball states cannot be distinguished in this sector The 
most interesting states to study there would be the 1^+ glueballs, that cannot be made of two gluons. 



D. Comments on gluelumps 

Gluclumps have been proposed in lattice QCD as a first approach to model gluino-gluon bound states [3l|. They 
are pure gauge states, as glueballs, but this time a static, scalar, color-octet source is present in the system. It appears 
that the gluelump masses are particularly low (the static source is decoupled from the system), actually lower than 
the lightest glueballs; the Hghtest gluelump is indeed a 1"' state of mass 0.87±0.15 GeV [32]- This can be understood 
in terms of constituent gluons because the presence of the color-octet source allows for only one gluon to be bound 
in a color singlet. Gluclumps seen as one-gluon states have been studied in detail within the Coulomb gauge QCD 
approach [33l |. but also with a spinless Salpeter Hamiltonian [33 |. The gluelump helicity states can be found in this 
last reference and read 

|7+;(J>1)^) = |0,1; J^,M, l) with P = (-)"', 

|T-;(J>1)^) =|0,1;J^,M,-1) with P = {-Y+K (7) 



They have a negative charge conjugation and are such that y^-^y = J{J + 1); formula ^ also holds in this case. We 

point out that using longitudinal gluons would lead to a low-lying gluclump, in disagreement with the lattice 
data. 

Interestingly, one-gluon states have been studied in lattice QCD more than twenty years ago in order to estimate 
the constituent gluon mass [35|. The following idea has been applied: The potential between two static color-octet 
sources should saturate at the energy at which a pair of gluons can be created from the vacuum. The energy of a 
confined gluon is thus half that energy. One arrives at an estimate of 500-800 MeV [HI. In the bag model, the same 
quantity can be computed and is given by 740± 100 MeV [36| , a value compatible with the 1"' gluelump. Interpreting 
the mass of that gluelump as a constituent gluon mass, one estimates the low-lying two- and three-gluon glucballs to 
be around 1.74 ± 0.30 GeV and 2.61 ± 0.45 GeV, in agreement with lattice computations 



III. THE RELATIVISTIC FLUX TUBE MODEL 



Once glueball helicity states are known, a Hamiltonian remains to be specified. To this aim, we present in this 
section the relativistic flux tube model. Although its original formulation, concerning mesons, relies on purely phe- 
nomenological grounds [13, l37l | , the relativistic flux tube model has been shown to be a well-defined approximation 
of the nonpcrturbative QCD interactions between a quark and an antiquark [isLlssll. We recall here the main aspects 
of this model, and refer the interested reader to these last two references for more details. 



A. Two-body case 



Let us begin by the case of a meson. The study of the quark- antiquark Green function in QCD involves the Wilson 
loop whose study eventually leads to crucial information about nonpcrturbative effects in QCD. In particular, an 
effective quark-antiquark Lagrangian can be obtained from the Wilson loop provided that the quenched approximation 
is made and that only the long-range, spin-independent, part of the interactions is considered. The spin- and short- 
range effects are actually not expected to be dominant with respect to the long-range interactions corresponding to 
the confinement. This covariant effective Lagrangian reads ^15. ,38j 



C = —mi 1 / x^ — 7712 V ^2 ^ a / dd \/ (w w')^ — w^ w'^, (8) 

with the boundary conditions 

w'^|,.i=x^ wn,=o=4- (9) 

The first two terms of Lagrangian ([8]) are kinetic terms corresponding to the quark and the antiquark seen as massive 
spinless color sources, whose position in Minkowski spacetime is given by the four- vector Xi(T). r is a timelike evolution 
parameter of the system, and the dotted quantities denote a derivation with respect to this parameter. The third 
term of Lagrangian ([S]) is a Nambu-Goto Lagrangian. It describes a string of generic energy density a linking the 
constituent particles through the boundary conditions We notice that the string coordinates w(0, r) also depend 
on a spacelike parameter labeling its different points; the prime denotes a derivation with respect to d. This string is 
the so-called relativistic fiux tube (or QCD string). It is a dynamical object arising from the nonpcrturbative part of 
the strong interaction. 

As mentioned in the introduction, the formation of flux tubes in QCD is supported in the meson and baryon cases by 
lattice QCD simulations. Moreover, it is of particular interest for the present study to stress that recent lattice QCD 
computations strongly favor the flux-tube picture in systems with constituent gluons like quark-antiquark-gluon (39j 
and three-gluon systems [i^]. Notice that, in all these lattice studies, quarks and gluons are seen as static color 
sources. Finally, Lagrangian ([8|) has been shown to be valid also for two-gluon systems in Ref. [ll|. The relativistic 
flux tube picture thus appears to be an interesting effective approach in the case of hadronic bound states. A comment 
on the value of a should be done at this stage. Let us denote a the fundamental string tension, i.e. the energy density 
of a flux tube linking a quark and an antiquark. Then, several lattice studies suggest that the adjoint string tension, 
that is energy density of a flux tube linking two static color-octet sources, should be given by ctq = (9/4)(t (bqI. l40l. l4ll| . 
In other words, we consider that the string tension scales as the quadratic color Casimir operator: This is the Casimir 
scaling hypothesis. 

It is convenient for further calculations to work in the temporal (or instantaneous) gauge, that is x" = X2 = r. In 
that gauge, the solution of the string equations of motion is simply given by the straight line w = xi + (1 — ^?) X2; 



see for example Notice that the hnearity of the flux tubes is in agreement with lattice QCD [l7|, |l8|. Once 

Lagrangian ([S]) is rewritten in the temporal gauge and with a straight flux tube, the relativistic flux tube equations 
can be obtained in the Hamiltonian formalism [isl . Is^ . In their quantized version, those equations read, in the 
rest frame and in the equal mass case (mi = m2 = m) \l4 . [2l| 



2 t ' V 

where {A, B} — AB + BA denotes the anticommutator, and where 



L^j = {vj,Wr} + a{rJ{v)}, (10) 

a ( arcsin w 1 , ^ 

H = {l,Wr} + -\r, , (11) 



Wr = ^M^\ ^ , ^^7t=^- (''^ 

V represents the modulus of the transverse velocity of the particles, and commutes neither with the radius r nor 
with the radial momentum pr- That is why various symmetrizations have to be performed. The complexity of these 
last equations comes from the fact that the contribution of the flux tube does not reduce to a static potential: This 
object is dynamical, and carries orbital angular momentum as well as energy. In consequence, v cannot be extracted 
from (fTUj) and the Hamiltonian cannot be written in a closed form. One can observe in (fTTj) that the linearly rising 
potential commonly used to describe confinement is supplemented by terms involving u, thus by relativistic effects due 
to the flux tube. Also the kinetic part is logically affected by these relativistic effects. We mention for completeness 
that equations generalizing (jlOp and (fTT ) can be obtained in the meson case when the spin of the quark and of 
the antiquark is no longer neglected [43|, |4J] . The resulting equations are quite complicated and we will not write 
them here explicitly. One actually gets couplings between the quark spin and the flux tube that can be understood 
as a Thomas precession in the color field generated by the flux tubes [i^. We point out that, for computational 
convenience, the following assumption will be made 



'L^)^^{L^). (13) 
Various limits of the equations (fTO|) and (|lip can be considered. When (i^) = 0, then one has trivially v = 



and the Hamiltonian of the problem simply becomes H — 2\/jP+rn? + ar since = + /r^ , that is a widely 
used spinless Salpcter Hamiltonian with linear confinement. In the nonrelativistic limit, an expansion in 1/m can be 
performed, and the resulting Hamiltonian is 2m + p^/m + ar, a Schrodinger equation with linear confinement. On the 

contrary, if m = and if {^^"^^ is defined as £{£ + 1), it can be shown by a semiclassical analysis that ~ 2TTai at 

large £ [S^j- These are the well-known Regge trajectories, observed in particular in light meson spectroscopy. Finall y, if 
the dynamical flux tube contribution is neglected in (fTO|) . one is led to the standard spinless Salpcter Hamiltonian [l4l | 



Ho = 2\/p^ + m2 + ar. (14) 

The corrective term 



a (L 



AH ft = TTT — r, with 



fir (6/i + ar) ' 



+ (15) 



can also be added in perturbation to take into account the dynamical contribution of the fiux tube at dominant order. 
The Hamiltonian Ho + AH ft defines the so-called perturbative flux tube model, that has already been studied and 
successfully applied to describe light mesons in (2^. |45|| . 

For two-gluon glueballs, one has obviously m — nig with nig the giuon mass. Moreover, we set a = Ua = (9/4)cr 

and ( i ^ ) is given by ^ . 



B. Many-body generalization 



The relativistic flux tube model can be generalized to many-body systems as follows. From the Casimir scaling 
hypothesis, one considers that each color source generates a straight fiux tube whose energy density is proportional 



to its quadratic color Casimir operator. Then, the flux tubes have to meet in one or several points such that the total 
energy contained in those flux tubes is minimal. In baryons for example, the junction point is called the Toricelli 
(or Fermat or Steiner) point of the triangle made by the quarks, leading to a Y-junction in agreement with lattice 
computations [H, EE S^l- For multiquark states (tetraquarks, pentaquarks,. . .), lattice studies indicate that the 
confining potential is compatible with a multi Y-junction flux-tube configuration [S, ^4911, again minimizing the total 
flux-tube length. The flux tubes then draw the so-called Steiner tree linking the quarks and antiquarks of the system. 
There may be more than one Steiner point in such cases: Numerical algorithms exist to compute them; see [50| for 
example. 

When gluons are present, the situation is different. Indeed, a constituent gluon may either generate an adjoint flux 
tube or two fundamental ones. In the two-gluon case, lattice results favor an adjoint string linking both gluons [39l.l4ll| 
rather than two superimposed fundamental strings. In a three-gluon system however, recent lattice results indicate 
that a triangle conflguration in which each gluon is linked to the two others by a fundamental flux tube is preferred 
to a Y-junction made of adjoint strings fioj. This result can be understood as follows from energetic considerations. 
Let us first consider that each gluon in a three-gluon system, located at the apices of a triangle ABC, generates 
an adjoint flux tube, and that the flux tubes meet at a point T. Assuming the Casimir scaling hypothesis, the 
energy contained in the Y-junction is then (9/4)(t(AT -|- BT + CT). But the triangular inequality leads us to 
(9/4)cr(Ar + BT + CT) > 2a{AT + BT + CT) > a{AB + BC + CA) [U: The energy of a triangular conflguration 
is energetically more favorable than the Y-junction, as also suggested by the lattice study [40l |. 

The energetic argument can be extended to a system made of Af gluons. Let us assume that the most favored 
configuration of a {Af — l)-gluon system is the one in which each gluon is linked to its nearest neighbors by a 
fundamental flux tube. We know that it is the case for A/" — 1 = 3 gluons. Then we choose a pair of gluons, denoted 
as X and Y, that each generate an adjoint flux tube. These flux tubes meet in Z, from which a flux tube starts and 
connect in some way to the Af —2 remaining gluons. Z is allowed to be in the 8, 10, 10, or 27 color representations. 
The less energetic flux tube generated by Z is obtained when Z is in a color octet. Consequently, Z and the Af — 2 
remaining gluons are equivalent to a {Af — l)-gluon system, for which the fiux tube configuration is known (closed- 
chain-like configuration). But at the same time, XYZ is equivalent to a three-gluon system: X, Y, and Z must be 
connected by fundamental flux tubes. One finally concludes that the most favored conflguration is the one in which 
Z is absent and in which a given gluon is linked to its nearest neighbors by two fundamental flux tubes. A schematic 
illustration of the above discussion is shown in Fig. [TJ 




FIG. 1: Graphical illustration of the most energetically favorable A/'-gluon configuration. 



This suggests that a relevant A/'-gluon Lagrangian generalizing ^ for A/" > 2 is 

f\e 7(w,w^)2-wfwf, (16) 
with 

Wi — 9xi + {1 - 9)xi+i and x^r+i^xi. (17) 

The resulting Hamiltonian equations are quite complex, see for example jisi] in the case of a Y-junction. Let us work 
again in the instantaneous gauge, and then neglect the contribution of the flux tubes to the momentum of the system. 
With this approximation one is led to the spinless Salpeter Hamiltonian 

A/' AT 
i=l i=l 

It is worth saying that not only this last Hamiltonian has to be considered in order to make explicit computations, but 
also the nontrivial symmetries of the wave function induced in the various J^*^ channels by the color wave function 
(commanding the charge conjugation) and the Pauli principle. These symmetries can become nontrivial as soon as 

A^> 3 m. 
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FIG. 2: Squared masses in units of the string tension a versus I for n = 0. A comparison is made between the results of 
the relativistic flux tube model (circles), the perturbative flux tube model (triangles), and the spinless Salpeter Hamiltonian 
(squares). The linear fits (|19p - (|21|l are also plotted. Calculations are made in the ultrarelativistic limit = 0. 



IV. REGGE TRAJECTORIES 



A. Perturbative flux tube and spinless Salpeter Hamiltonian 



So far we have presented three different types of Hamiltonians that can be used in hadronic physics. The first one 
is the relativistic flux tube Hamiltonian pH)) . that contains the full relativistic dynamics encoded in Lagrangian ([S]), 
itself derived from QCD. Numerical computations become simpler if one is able to deal with a closed Hamiltonian; 
that is why the perturbative flux tube model, defined by (|14p and (jl5p . can be used as an approximation of the full 
relativistic flux tube model. Finally, the dynamical contribution of the flux tube may completely be neglected and 
the resulting Hamiltonian is the spinless Salpeter one ([H]). 

In this section, we compare the three frameworks in the ultrarelativistic limit, i.e. mg = 0, that is the case in 
which the relativistic effects are the most important. It is well-known that in the case of massless particles, all those 

approaches lead to particular orbital (Regge) and radial trajectories: Defining {^"^^ — ^(^ + l)i it is observed that 

(X i and (x n for large ^ and n respectively. Notice that, in the meson case, i is the orbital angular momentum 
of the system, while i{l + 1) = J( J + 1) + 2Ai A2 for the glueballs. Let us now give a more quantitative description of 
that property. 

The only energy scale in the massless case is \fa. Consequently, hP' ja is a dimensionless quantity that can be 
computed numerically. Our procedure is the following. Thanks to the Lagrange mesh method, which is a simple 
and accurate numerical technique, we are able to solve the relativistic flux tube model [5^ and spinless Salpeter-type 
Hamiltonians [ssj . Then, we compute the square masses for each approach and for < €,n < 8, which is enough 
to clearly observe the expected linear trajectories. As an example, we plot in Fig. [2] the n = Regge trajectory 
computed in the three considered approaches. The linear behavior is clearly observed, but the Regge slope is different 
in each case. Fits on the complete sets of data that we have computed show that the square masses are in excellent 
agreement with the following linear forms: 

= 27ra(1.944n+ 1.238£+ 1.606) (19) 

for the spinless Salpeter Hamiltonian (fT4| . 

Af2 = 27ra(1.968n+ 1.092 £+ 1.586) (20) 

for the perturbative flux tube model (fTi|) - p3|) . and 

= 27ra(2.003n+1.018£+ 1.533)«27ra(2?i + £ + 3/2) (21) 

for the relativistic flux tube model (llO |l - (fTT|) . The qualitative behavior of the mass spectrmn in all those three 
approaches is clearly identical. That is why, even if a spinless Salpeter Hamiltonian is a crude approximation of the 



relativistic flux tube model, it can reproduce the same data provided that the parameters are correctly fitted. Indeed, 
let us imagine that we want to reproduce the Regge trajectories of light mesons. Since the Regge slopes (in units of 
a) are quite different, 7.779, 6.861, and 6.396 for equations (|19 p - (|2T|) . one has to choose a different value for the string 
tension in each case in order to agree with experiment. 

It is also worth mentioning that the radial and Regge slopes are never equal, even approximately. It turns out from 
the fit ()2ip that the relativistic flux tube model actually exhibits a remarkable harmonic oscillator degeneracy, with 
a radial slope twice as large as the Regge slope. It implies in particular that, when applied to light mesons, the flux 
tube model alone is unable to explain the Coulomb-like degeneracy that seems to be experimentally observed in the 



highly excited mass spectrum [5J|. Notice that the slope of radial trajectories is almost insensitive to the various 



approximations, while the Regge slope varies significantly. This is in agreement with 55], where it is shown that the 
relativistic flux tube model is equivalent to a spinless Salpeter Hamiltonian for large radial excitations. 



B. Open-closed flux tube duality 

It is readily checked that the mass spectrum (f?T|) corresponds not only to the relativistic flux tube model but also 
to the following Hamiltonian 

He = + (27rar)2 , with ^ = ^- (22) 

This last Hamiltonian has a nonstandard form, which actually corresponds to a closed circular flux tube of energy 
density a = a/2 and of radius r l/\/5 [HHl- Consequently, the following models lead to equivalent mass spectra: 
An open flux tube of energy density a and length r with massless particles at his ends and a closed circular flux tube 
of energy density a/2 and radius r. 

Such a duality is particularly interesting in the case of glueballs since both the closed- and open-flux tube approaches 
can be found in the literature. The open flux tube framework has the advantage of providing a common description of 
light mesons and glueballs [2^. But, it is sometimes argued that gluons in a glueball should all have the same status 
and not being either constituent or part of a flux tube. That is why some authors have proposed that glueballs should 
be described by closed fundamental flux tubes 0, [H, [57| . Such a picture has also the advantage of being closer to 
the lattice QCD formulation, in which pure glue states are described by loops on the lattice [53] ■ In the closed flux 
tube approach, the different J^^ states are generated by orbital or phononic excitations of the flux tube. It is worth 
noting that, if we set a — aa for the open flux tube in a glueball, a/2 is not far from the expected value a. 

For completeness we can quote that more exotic, knotted, configurations of the flux tubes have been suggested 
to be relevant in the modelization of glueballs [HI], and that a description of glueballs as two-dimensional spherical 
membranes has also been proposed in [5^. In this last reference, the effective Hamiltonian reads \/p'^ + {jiirTr'^Y, 
the potential term being the surface of the spherical membrane of tension T . The squared mass o per ator is in this last 
case a quartic oscillator, whose spectrum would be of the form M^/T^/'^ (x (2.27n-|-^-f 1.68)^/'^ [60|, in disagreement 
with the open flux tube behavior. 

The A/'-gluon generalization (fTB|) of the flux tube model is nothing else than a closed fundamental flux tube with 
transverse constituent gluons attached on it. That picture seems roughly compatible with an open- closed-flux tube 
duality, with the dynamics of the constituent gluon seen as an alternative way of modeling the various phononic or 
orbital excitations of a pure closed flux tube. However, we point out that the celebrated closed flux tube model 
of Isgur and Paton [57| is not in good agreement with recent lattice data when taken in its original form [gj]. In 
particular, it has difficulties to reproduce the important splittings between C — ~\- and — glueballs with the same . 
It thus appears that improvement of this model are needed in order to understand the lattice QCD glueball spectrum. 



V. THE GLUEBALL SPECTRUM 



The relativistic flux tube model presented in section IIIII provides us with a relativistic - although not explicitly 
covariant - description of hadronic bound states. Supplemented by other physical mechanisms that we present 
hereafter, it can be used to understand the lattice data. 



A. Short range interactions 



The flux tube stands for the long-range interactions, i.e. the confinement. At short distances however, other 
physical mechanisms become dominant. These are mainly given by the one-gluon-exchange processes between color 



sources. The corresponding short-range potentials can be computed from the QCD Feynman diagrams at tree level or 
from the perturbative part of the Wilson loop. In the quark-antiquark case, one is led to the well-known Fermi-Breit 
interaction. Other cases such as gluon-gluon and quark-gluon potentials have been investigated for example in [63 |. 
As expected, it appears that, at the dominant order, the short-range potential is always given by a Coulomb term 
scaled by some color factor depending on the considered color sources; in particular one is led to 

^Hoge = -3— (23) 
r 

for a two-gluon glueball, as being the strong coupling constant and r the distance between the two gluons. In a 
three-gluon system, the 3 factor is replaced by 3/2 because a given gluon pair is in a color octet. In the special case 
of the four-gluon 0"' state, gluon pairs are found to be either in a color decuplet or antidecuplet, leading to the 
vanishing of AHoge [111 . 

Already at one-loop, as is running with the energy scale q^; we know that a5(q^ — > oo) = 0. Moreover, we consider 
that a5(q^ — > 0) = ao. This is coherent with lattice QCD calculations, showing that the static potential between a 
quark and an antiquark is accurately fitted by the form ar — (4/3)as/r [63l |. An intuitive justification of this last 
relation is that the total energy of the system is separated into a dominant, nonperturbative, part (the flux tube), 
and a residual, perturbative. part (the o ne-g luon-exchange potential). This is in agreement with the aforementioned 
lattice works [ll, [l^, [H, Ho, lH, \^ [H, |43] , in which the potential energy between the various sources considered 
is always compatible with the expected one-gluon-exchange term at short distances, while the flux-tube contribution 
comes into play at long distances. 



B. Instanton-induced forces 

The influence of instantons on light quarks and light mesons has been intensively studied (see ^64| for a complete 
review about instantons in QCD). It is known, in particular, that instanton-induced forces in the lightest mesons can 
be included in potential quark models as an isospin-dependent term whose contribution is nonzero in the pseudoscalar 
channel only [65j . Such a term is needed to reproduce the very low value of the pion mass without lowering the masses 
of the other mesons. 

Unlike the meson case, the effects of instantons on glueballs have been less intensively studied. In a pioneering 
work [g^ , it has been computed that direct instantons induce a strong attractive force in the scalar glueball channel 
and a strong repulsive force in the pseudoscalar one, while no effect is observed in the tensor channel and presumably 
in the other channels. It is moreover tempting to deduce from (67| that these forces in the scalar and pseudoscalar 
channels are of equal magnitude but of opposite sign. This has led us to propose the following phenomenological 
instanton mass term in Ref . [l^ : 

AHi = -PISj^o, (24) 

where P is the parity operator and where the arbitrary parameter X is assumed to be positive and constant in first 
approximation. It is relevant to include instanton-induced forces in glueballs because they play an important role in 
the meson case and, since they are needed in glueballs as well jB^, i67|], it is coherent to include them in the model. 



C. Two-gluon spectrum 



Using the Lagrange mesh method , we are now able to compute the two-gluon glueball mass spectrum from (fTU|) 



and (fTT|) in which H is replaced hy H + AHoge + AHj and in which (L^) is computed thanks to ©, coming from 
the helicity formalism. The corresponding J^'~'' states have all a positive charge conjugation and are given by ([5|). 

The parameters remain to be fixed. We set nig = in the spirit of [l^ . which was a first attempt to include the 
transversality of the constituent gluons in apotential model. We think that it is interesting for the present model 
parameters not to be in disagreement with [22] , in which the experimental mass spectra of light an heavy mesons are 
quite satisfactorily reproduced within the perturbative flux tube model. We thus keep the value as = 0.4, used in this 
last reference in the light meson sector. The value used for the string tension was 0.185 GeV^ in Ref. [22], ensuring 
a correct reproduction of the experimental mesonic Regge slopes. Let us denote tXp/t = 0.185 GeV^ the fundamental 
string tension fitted on light mesons in the perturbative flux tube model jUj. Then the same Regge slope will be 
obtained by using a = (6.863/6.398)crp/t — 0.198 GeV^ in the relativistic flux tube model. In summary we can set, 
in the two-gluon relativistic flux tube model 



mg = 0, aa = (9/4)(T, a = 0.2 GeV^ as = 0.4, and I = 0.450 GeV, 



(25) 
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TABLE L Comparison between the C = + glueball masses computed in pure gauge lattice QCD [2,[3,[63| and the spectrum of 
the relativistic flux tube model with parameters (|25p . The corresponding helicity states are given in the last column. 



where the magnitude of instanton-induced forces is fitted on the lattice splitting between the scalar and pseudoscalar 
states. 

The mass spectrum we obtain by computing two-gluon glueball masses within the relativistic flux tube model is 
given in Table [J It can be seen that the agreement with the states that are currently known in lattice QCD is good, 
all the computed masses falling into the error bars. Only the 3"'""'" state is missed; we indeed find a very low value 
for it, so low that the lattice mass rather corresponds to the first radially excited state. To our knowledge, there is 
no theoretical reason to exclude the S"'"^ state with n = 0. Since that glueball is the only one coming from the £>_ 
family, it would be interesting to compute the 5"^'^ glueball mass on the lattice and see whether it agrees more with 
the n = or n = 1 state in our model. 

It is worth pointing out that the masses of the states (J > 0)+) and (J > 0)^) arc equal in our approach 
because of ([T3| . It was already suggested in the pioneering work '8] that such a characteristic pattern would be an 
interesting tool to identify glueball resonances. Moreover, this degeneracy cannot be lifted by the introduction of 
tree-level spin-dependent operator like spin-orbit, spin-spin, ttc. operators [sl. [l^. A mixing term between the 5+- 
and i'^-states would be able to affect this mass degeneracy. The existence of such a mixing can be understood as 
follows. Let us consider two transverse gluons on mass shell, with their helicity denoted by + or — . If ^(Ai, A2; A'j^, A2) 
is the gluon-gluon scattering amplitude with A^ and A^ the incoming and outcoming helicities respectively, then it can 

be shown that, at tree level [fiil, yi(++; ++) ^ and ^(++; ) ^ 0, but ^(++; H — ) = — V) = 0. In other 

words, the helicity-conserving amplitudes are nonzero at tree- level and give rise at the dominant order to the Coulomb 
term ([23]). However, the scattering amplitudes of two gluons is zero when only one helicity is reversed. Reversing 
a single helicity transforms a helicity-singlet in a helicity-doublet and vice-versa. It is thus logical that we obtain 
no singlet-doublet mixing with tree-level relativistic corrections. At the next order however, scattering amplitudes 
such as A{-\ — h; H — ) become nonzero and a singlet-doublet mixing can appear. However, the explicit form of the total 
scattering amplitude beyond tree level is very complex (see Ref. [69l|). and computing an effective potential from all 
the allowed diagrams is a task which is beyond the scope of this paper. Another consequence of such a mixing is that 
it decreases the mass of the lowest-lying {2k + 2)++ states, improving in particular the agreement of the 2++ state 
with the optimal lattice value. 

The glueball spectrum of Table U has been computed by considering massless constituent gluons. It is worth 
mentioning that other frameworks like Coulomb gauge QCD rather use a nonzero gluon mass of typically 600 MeV [|5]. 
We have checked that a mass spectrum in agreement with lattice QCD can be obtained with nig = 0.6 GeV, a = 
0.135 GeV^, as = 0.45, and I = 0.45 GeV, but the global agreement is less satisfactory than with parameters (|25p . 
Those values are inspired from a Coulomb gauge study of hybrid mesons (70j . 

Finally, it is worth mentioning that the auxiliary field method allows to find an approximate analytical mass formula 

for the Hamiltonian 2y^^^2~fm^ + ~ 3as/r — PXSjq, which is a simplification of the full relativistic flux tube 
model [7l[. Correcting the formula obtained by taking into account the result (dH), the final mass formula, for 



TOg < (Ta, reads [nl 



where 



M2g = Mo + —rT^ P^Sj^o, (26) 



Mo = ^yTraa{4n + 2e + 3{l-as)) (27) 



and where n and £ are the radial and orbital quantum numbers of the gluon pair. £ is computed through the relation 
£{£ + 1) = J( J + 1) + 2A1A2 [see Eq. This last formula reproduces the glueball spectrum for two massless gluons 
computed within the relativistic flux tube model with a global accuracy better than 4% (but 9% at J = 0). 



D. Three- and four-gluon systems 



It is possible to obtain an estimation of the lightest three-gluon glueball mass thanks to analytical results obtained 
with the auxiliary field method in [73 |. It is shown in this last reference that an accurate approximate analytical 
mass formula for the Hamiltonian H^g = + X]?<j=i i^l^i ~ \ ~ (3/2)a5/|a;i — intended to describe 

glueballs made of three gluons, is given by 

2 

Mig ^ {I2V3 N - 5Aas)cr, with N = 3 + ^{2ni + £i), (28) 

2 = 1 

where Ui and £i denote the radial and orbital quantum numbers associated to the two Jacobi coordinates. This last 
formula is actually an (accurate) upper bound of the exact mass spectrum [72]. We also point out that it has been 
obtained within a spin-independent formalism. In other words, only the orbital angular momenta are present, and no 
input coming from the gluons hehcity has been added. The ground state of (|28l) has a mass equal to (36-\/3 — 54as)cr, 
that we expect to be an underestimation of the ground-state mass. Indeed, that value is obtained by setting all the 
relative momenta equal to zero, i.e. N = 0, and it can be seen, already in two-gluon case, that the introduction of 
transverse degrees of freedom adds components of higher orbital momenta in the system, tending to increase the total 
mass [13] ■ With the parameters we obtain Al^^g = 2.855 GeV; since the lightest C = — glueball computed on the 
lattice has a mass of 2.980 ± 170 GeV, it appears to be a relevant lower bound. 

We proposed in Ref. j23] to see the 0"' glueball as a low-lying four-gluon state rather than a highly-excited three- 
gluon state. Let us check that hypothesis. The particular color structure of this state (see section III C|) causes the 
short-range interactions to vanish, and enforces a particular mixed symmetry that suggests an approximate way of 
computing the 0^ mass by using two-gluon clusters (see Ref. [l^] for more details). First, we compute the mass of the 
lightest antisymmetric state made of two gluons linked by a fundamental flux tube of tension cr. This state is a J = 1 

one, for which {^^^^ = 4 following and the relativistic flux tube model gives a cluster mass of Mc = 1.964 GeV. 

Second, we couple both massive clusters in a S'-wave channel, which is the lightest symmetric configuration of these 
clusters, linked by a flux tube of energy density 2a. The final mass is 4.918 GeV, a value compatible with the lattice 
one of 4.780 ± 0.290 GeV Q . 



VI. GLUEBALL DECAY WIDTHS 



Together with the mass, the total decay width is a particularly relevant observable characterising a given hadron. 
We adapt in this section the method of [73\, giving a way to estimate the two-gluon glueball decay widths within a 
flux tube approach. 



A. General mechanism 



Let us consider a two-body hadron (meson or glueball typically). The flux tube linking both constituents is a 
particular configuration of the gluonic field in QCD; it is purely chromoelectric in the rest frame. Within this 
chromoelectric field, a particle-antiparticle pair can be created from the vacuum thanks to the Schwinger effect [74| . 
This creation eventually leads to a breaking of the flux tube, and thus to the decay of the considered system. In 



the case of a meson, say QQ, a quark- antiquark pair qq can be created inside the flux tube and causes the meson 
to decay as QQ — > qQ + qQ. It is known for a long time that such a mechanism can successfully describe the decay 
widths of mesons, and even baryons [t^. The idea can be generalized to the glueball case, with the creation of a 
quark-antiquark pair within the adjoint flux tube. Other processes can also occur in glueballs: Simultaneous creation 
of two quark-antiquark pairs within the flux tube, quark-antiquark disintegration of a constituent gluon, etc. Many 
processes are actually allowed because more energy is stored in an adjoint flux tube than in a fundamental one, and 
glueball decays are clearly more complex than meson ones. That is why the calculations we will perform here, which 
are based on the qq creation picture, only aim to give rough estimations of the glueball total decay widths. We remark 
that this approach should not be used in the case of heavy quarkonia since decays induced by short-range processes 
may become dominant with respect to the flux-tube induced ones. 

The quantity (r^) is a measure of the flux-tube length. It is plausible that the probability to produce a quark- 
antiquark pair is proportional to the string length since this production is linked to the energy carried by the chro- 
moelectric field. The energy scale of the flux tube model with string tension a is given by in the ultrarelativistic 
limit. We are thus led, in agreement with [73|], to the following ansatz for the decay width 

T = jV^i^/{^-xo), (29) 

where x = r^fa is used together with q = pj \fa to get dimensionless average values, and where 7 and xq are assumed 
to be universal constants. Notice that (|29|) could be written on the alternative form F — 7(M — xq), where the hadron 
mass explicitly appears. 
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FIG. 3: Plot of the spin- and isospin-averaged experimental decay widths of nn mesons sharing the same spectroscopic assign- 
ment (n + \)^^^^Lj versus ^ (x^) (circles). The form l|29|l with the parameters (|30p and a = a is plotted for comparison (solid 
line) . 



B. Numerical results 



One has first to fit the parameters 7 and xq. To do that, we resort to the available experimental data concerning 
light mesons [TG^, that we denote as nn mesons, n standing for u and d quarks. The present approach is spin- and 
isospin-independent. Consequently, we have to apply it to the spin- and isospin-averaged decay widths of the lightest 
mesons, as done in [t^. Only the lightest states are taken into account because they are expected to possess the 
largest radii, thus to be the ones for which flux-tube induced effects are the most relevant. The values ^ (x^) are 
computed by using the relativistic flux tube model for the quark-antiquark pair with the Coulomb term — (4/3)as/a; 
and the parameters ([^5]) . The couples (a/ (x^), F) are then fitted to (|29p with a = a, leading to 

7 = 0.300 ± 0.045, xo = 0.827 ±0.309. (30) 

The result in the light meson sector is plotted in Fig. [31 
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TABLE II: Total decay width of some low-lying glueballs computed with formula (|29|l . the parameters (|30p . and a = da. The 
error bars on F^j come from the uncertainty on 7 and xq. The corresponding masses are recalled in the last column (see 
Table m. 



We can then estimate the total decay widths of some glueballs by computing the corresponding values of {x^) ■ 
Some results are shown in Table |lll The decay widths of the scalar and pseudoscalar glueballs are degenerate in 
our approach, since instanton-induced interactions are only added as a constant term and do not affect the radius. 
It is interesting to compare these results with some experimental data. It has been suggested in several works 
that the /o(1710) could be a mainly glueball state [t^]- Its experimental mass and decay width are respectively 
1.724 ± 0.007 GeV and 137 ± 8 MeV The decay width is lower than our estimate, although still falling in the 
error bars of our prediction. It is worth mentioning the detailed Coulomb gauge computation of 78] leading to the 
quite low value of 100 MeV for the scalar glueball decay width. Those results could favor the glueball interpretation of 
the /o(1710), although this last work was rather concerned with the /o(1810) resonance. There are two tensor states 
which are located in the mass range of lattice QCD: The /2 (2300) with a mass and decay width of 2.297±0.028 GeV and 
149 ±40 MeV, and the /2(2340) with a mass and decay width of 2.339 ±0.060 GeV and 319t?o MeV respectively 
The masses of these states are compatible, but their decay widths are not. Our results rather favor the /2(2340) as a 
tensor glueball candidate. 



VII. THE PSEUDOSCALAR GLUEBALL 



Recently, the "q-r]' mixing angle as well as the gluonium content of the 77' have been computed by the KLOE 
collaboration from a study of the channels 0(1020) 777 and r]'^ |79j. Let us recall the results obtained in these 
last references. First of all, the rj and t]' wave functions are are decomposed into three Fock-space components: 
\nn) = {\uu) ± \dd))/\/2, \ss) and a glueball component that we denote \gg). With the following parameterization, 

I?]) = cos (^p I nn) — sin I ss) , 
I77') — co&ipgSviiifp\nn) + cosipg cosipp\ss) + 

+ sin^g|55), (31) 

one obtains from a fit to the experimental data [t^ 

<Pp (40.0 ±0.7)°, Zl ^ sinVg = 0.13 ±0.04 

^ (ySg = (21.1 ±3.4)°. (32) 

It is assumed that the 77- resonance has no glueball component, while the gluonic content of the r/ is found to be 
nonzero. 

The parameterization ([31]l allows the existence of a third state, orthogonal to the ry and 77' quantum states, expressed 

as 

\G) = — shiLpg siwipp \nn) — sinipg cosipp \ss) + cosipg \gg) . (33) 
Introducing rotation matrices, it is readily checked that 

\r])\ /l ^ \ /cos(/3p —siTupp 0\ f\nn)\ /|rin)\ 

W) \ ^ \ ^ cos (fig simpg sinipp cos (pp \ss) = U \ss) . (34) 
JG)/ \0 -sin(^g cos^^y V 1/ \\gg) J \\gg) ) 



Let us now define M = diag(M^, Af^/, Mq)- Tlien, tlie states (IMl) are eigenstates of tlie mass matrix 

M = UUfU. (35) 

The diagonal elements of M are the masses of the unmixed pseudoscalar |nn), |ss), and \gg) states; we denote Mgg 
this last element. Mgg can be a priori different from the pure gauge pseudoscalar glueball mass since quark loops 
inside the correlator may change its value. A possible effect of sea quarks is to decrease glueball masses (see for 
example Ref. [13!] for more information). In the scalar glueball mass with Nf = 2 is found to be 15% lower than 
the pure gauge mass. But no definitive conclusions can be drawn at this stage (see Ref. [8l|, in which non significant 
difference with the pure gauge mass is found). 

^From the results p2p . \G) appears as being a glueball at 87%; it could actually correspond to the "physical" 
pseudoscalar glueball. Its mass can be estimated from the KLOE data. Some algebra indeed shows that the following 
relation holds 

Mg = Mgg + (Mgg - M,,,) tiill^ <f g . (36) 

Formula p6p gives the mass of the physical pseudoscalar glueball in terms of the experimentally measured mixing 
angle ipg and the bare pseudoscalar glueball mass. Taking into account the error on the determination of ipg [79l |. we 
find 



Mg = (1.149 ± 0.056)Mgg ~ (0.143 ± 0.053) GeV, (37) 

where M^/ = 957.66 ± 0.24 MeV Using our pure gauge value Mgg = 2.610 GeV, we are eventually led to suggest 
the existence of a pseudoscalar glueball candidate with a mass of Mg = 2.856 ± 0.199 GeV and a total decay width 
around 220-315 MeV (we take the upper values of Table Ull because of the repulsive instanton interactions and because 
the mesonic components allow for more decay processes than in the case of a bare glueball). Interestingly, such a 
mass is not far from the T^c-meson, whose mass is 2980.3 ± 1.2 MeV but whose decay width is only 26.7 ± 3.0 MeV. 
For completeness, we mention that following a previous work [s^ . the Ty^- glue ball mixing could be large, while there 
are preliminary lattice results suggesting that the rj'-rjc mixing is small [83[. So far it is not known how large is 
the influence of sea quarks on the pseudoscalar glueball. In particular the 7y(1405) is sometimes considered as a 
pseudoscalar state with a large glueball component, hence it could correspond to \G). This conclusion is shared by 
the recent works [s^, IsHl , based on the FKS formalism 86] , with the same experimental input as in the present study. 
But following ([5^ . this would imply Mgg 1.35 GeV, a value that would demand huge sea quarks effects to be 
reached. 

It is finally worth saying that there is a controversy about the KLOE results [7§]. For example, one can quote 
Ref. |87[, in which it is found through another analysis of radiative scalar and pseudoscalar meson decays that 
(Pp = (41. 4± 1.3)° and (fg = (12± 13)°. Both the rj and rj' resonances are compatible with a zero gluonium component 
in this case. Notice however that more recent computations of the same author agree with a nonzero gluonic component 
in the ry' (ssl ]. 



VIII. GLUEBALLS AND THE POMERON 



It is well known that in the many high energy reactions with small momentum transfer, the exchange by the 
highest-lying Regge trajectory, called as a soft Pomeron, gives the dominant contribution. This exchange carries 
vacuum quantum numbers PC = H — h and has very peculiar properties in comparison with the usual Regge pole 
trajectories. From the analysis of pp and pp cross sections it follows that the Pomeron trajectory has a linear behavior 
reading [89|] 

a(t) = 1.08 + 0.25 i. (38) 

It is expected that such a trajectory should contain physical hadron states, with t = and ap{t) — J. But the 
Pomeron does not seem to be related to light mesons since the latter have lower intercepts and larger slopes. There 
has been a long-standing speculation that the physical particles on the Pomeron trajectory might be glueballs. In the 
lattice study [J], it is shown that the hghtest (2, 4, 6)++ pure gauge states lie on a Regge trajectory that is consistent 
with the Pomeron ([55]) . However, the masses that are found in Q are significantly lower than the lattice studies 0, Q 
on which we rely in the present approach. Furthermore, in a recent calculation, Meyer updated his value for the scalar 
and the tensor glueball masses [90|]. Although being obtained with a different method, the two new masses are, in 
this reference, in agreement with Refs. [1, Q]- We thus conclude that this last mass spectrum is a reliable one and 
that it is worth discussing a bit the compatibility of our results with the Pomeron. 



First of all, the linearity of the trajectory (IHS)) is in qualitative agreement with the flux-tube picture for glueballs. 
The lightest glueball states with J^'-^ =(even> 2)+"'" are generated by the family. Within our model, it is readily 
computed that the squared masses of this family, with J running from 2 to 20, can be fitted by the following form 

J = 0.34M2 - 0.41 « (Af2 - 1.16). (39) 

2TTaa 

The slope of this trajectory is quite different from 0.25 GeV~^; actually we should have taken the large value a = 
0.29 GeV^ to reproduce it. This is due to the relativistic effects contained in the flux tube. Within a spinless Salpeter 
formalism, for the same string tension, we would have found a slope of [7.77(9/4)0.2]"^ GeV~^ = 0.29 GeV~^, in 
better agreement with ([38| . One has thus to be careful when claiming that a potential glueball model reproduces 
the Pomeron slope: This is strongly model-dependent, and nonstandard values of the string tension can be needed 
to match the slopes. We also see in ((5^ that the intercept is not close to 1: Our glueball trajectory thus do not 
correspond to the leading Pomeron trajectory. It is worth mentioning that a large- J expansion of formula with 
nig — gives 

A/^ 3 

^ = ^ + 2 («s - 1) + O {M^g) « 0.35 M% - 0.9. (40) 

The slope is nearly equal to the fitted form, while the intercept is not. With that formula, J -I- 2 could agree with the 
Pomeron trajectory; this is a supplementary reason to consider that our glueball spectrum rather lies on a daugther 
trajectory. 

Let us now quote the work [9li] , where not only pp and pp cross sections are fitted within the Pomeron framework, 
but also TT^p and K^p cross sections. The following results are found: With a single (soft) Pomeron, fits are consistent 
with 



a(i) = 1.09 + 0.33 i. (41) 

With two Pomerons (soft and hard), fits are consistent with 

a^(t) = 1.07 -f 0.30 i and auif) = lAh + Q.lt, (42) 

although the need for a hard Pomeron is not crucial for the quality of the fits [9l|. These results lead to a soft 
Pomeron slope that is in better agreement with our slope, especially (|41]) . We can also mention that our 2++ mass is 
compatible with the point auC^), while in Ref. [l^l the tensor glueball is found to lie on the soft Pomeron trajectory. 

The fit ([55)1 would require the tensor glueball to have a mass around 1.92 GeV, while the fit (|¥T|l favors a tensor 
glueball mass around 1.65 GeV. Both are values far from our results. In conclusion, the conjecture that the lightest 
(even-J)++ glueballs lie on the Pomeron trajectory is not favored by our results; we would rather conclude that 
these glueballs are located on a daughter trajectory. But, we focused here in reproducing pure gauge lattice data, 
by computing the spectrum of genuine two-gluon states. It is plausible either that the leading Pomeron trajectory 
contains no pure glueball but rather mixed states (glueball, mesons, tetraquark,. . . ), or that the sea quark effects lead 
to a strong decrease of even the pure glue states. Both mechanisms escape our current understanding. 



IX. FINITE-TEMPERATURE EFFECTS 



The glueball lattice spectrum, given in Table |T] and satisfactorily reproduced by the relativistic flux tube model, 
was formally obtained at T = 0. It is well-known however that nontrivial phenomena occur at nonzero temperature. 
In particular, there exists a critical temperature in QCD, denoted as Tc, above which deconfinement takes place, 
eventually leading to the celebrated quark-gluon plasma. Since the relativistic flux tube model is an effective approach 
describing the conflnement in QCD, it should not be used at temperatures higher than Tc. However, it can be expected 
that, starting from T = and increasing the temperature, the flux tube will "weaken" due to thermal effects and then 
will disappear at T^. One can think of a progressive disorganization of the chromoelectric field because of random 
thermal fluctuations, as illustrated in Ref. [93|. That picture is supported by lattice computations of the static quark- 
antiquark energy, showing that the dominant effect below Tc is a diminution of the effective string tension, that can 
be parameterized by [93| 

a{T) = a[l-^y (43) 



The variation of the Coulomb term can be ignored in a first approach [93|. For completeness, we mention that a 
potential model describing glueballs above Tc can be found in Ref. [95^1. 

The behavior of the glueball spectrum at finite temperature has been first investigated on the lattice in Refs. [o^lsj]- 
The basic conclusion of those works is that the scalar glueball mass remains constant from T = to Tc up to 
fluctuations of order 100 MeV. Notice that such a behavior was expected from calculations relying on effective low- 
energy lagrangians [ol]. The decay width of the scalar glue ball, however, has been shown to increase significantly 
because of thermal effects [o^, [sl] . In the more recent work jlO0l | , the pseudoscalar and tensor glueball masses are 
also shown to be constant from T = to Tc, again up to an accuracy of 100 MeV. 

How can those results be qualitatively explained within the relativistic flux tube model? First, let us consider that 
the dominant contribution of the thermal effects is to modify the string tension according to Eq. (|43p . Second, the 
influence of T on the instanton-induced effects is neglected. Thus, in order to get a nearly constant scalar glueball 
mass, TOg has to increase with T in order to compensate th e dec rease of the string tension. The idea that gluons should 
gain a thermal mass has already been mentioned in Ref. [lOll | for example. Since no conclusion can be drawn about 
the form of mg{T), we consider it as a free parameter and compute it at various temperatures so that dTMQ++ = 
between T and Tc. The gluon masses numerically obtained are well fitted by the following form 



We now turn to the thermal decay width of the scalar glueball. We gave in Sec. I VII a simple way to estimate 
the glueball decay width in the case where the decay by pair creation is the dominant mechanism. Using the fitted 
form (jU]) and the formula it can be computed that \/Jr^ significantly increases from 2.850 GeV~^ at T = 
to 3.823 GeV~^ at T = 0.9 Tc. However, a(T) quickly decreases, so that a(T) y/ (r^) becomes negligible at high T, 
meaning that the decay by pair creation is no longer dominant. Wh at we have instead is a progressive "dissolution" 
of the glueball in the medium, expressed by the increase of (r^) with T; such a phenomenon should contribute 
to increase the decay width of the glueball. Moreover, since the constituent gluons become quite heavy near Tc, 
perturbative mechanisms such as those found in heavy quarkonia might come into play. It is known for exam ple that, 
because of in-medium effects, heavy quarkonia gain a thermal decay width that roughly grows like T jl02l |. Again, 
such a mechanism could increase the glueball decay width. We can thus conclude that an increase of the glueball 
decay width with T is not excluded a priori within our effective approach, but further theoretical investigations are 
needed to clarify and confirm that statement. 

To our knowledge, no widely accepted theoretical derivation of the thermal gluon mass between T and Tc can be 
found in the literature so far. We have shown that the form ([^^ leads to the expected features for the 0'*'"'' state. 
An interesting point to mention is the following: Only a state-dependent thermal mass is able to keep all the glueball 
spectrum constant with respect to the temperature. It can be seen from the approximate mass formula (|26p indeed 
that the whole glueball spectrum will have a constant mass provided that mg(T) is approximately given by 



Mq being given by Eq. (j27p . which moreover gives an a posteriori justification of the form (|44p . Further lattice 
studies concerning the evolution of the glueball masses versus the temperature in various channels could provide more 
information about the thermal gluon mass. If all the glueball masses appeared to be constant, then a state-dependent 
gluon mass like (|45|l would be needed. On the contrary, let us suppose that mg{T) is state-independent, but such 
that the scalar gluebal has a constant mass. Then, the pseudoscalar glueball will have a constant mass following our 
approach, since instanton-induced interactions are assumed to remain constant with the temperature, while the mass 
of the other glueball states will decrease with T. Indeed, the thermal mass is not large enough to balance the decrease 
of the string tension in that case. For example, using Eq. (|44l) . we compute that the 2++ mass could decrease from 
2.529 GeV at T = to 2.326 GeV at T = 0.9 T g- For the moment, lattice data do not allow to conclude that the tensor 
glueball mass significantly decreases; see [l00j |. Computations with increased accuracy might be quite interesting in 
order to check whether that behavior is observed or not. 




GeV. 



(44) 




(45) 



X. CONCLUSIONS 



The relativistic flux tube model has been firstly applied to compute masses and decay widths of glueballs seen as 
bound states of transverse constituent gluons. It allows not only to understand to qualitative features of the lattice 
spectrum, but leads also to a good quantitative agreement with these data. Figure U] is a summary plot of the masses 



we obtain. Notice a peculiar result we have found: An open flux tube with massless particles attached at his ends has 
the same mass spectrum than a closed flux tube with an energy density two times smaller. There exists thus a kind 
of duality between open- and closed- flux tube models, that draws a bridge between apparently distinct approaches 
existing in the literature. 

Only a few calculations of glueball decay widths exist nowadays (see Refs. [T^ItII and references therein), and are 
mostly concerned with the scalar glueball. Our framework allows to estimate the decay width in various channels, 
leading in particular to the conclusion that the /2(2340) resonance could be a relevant candidate for the tensor 
glueball. Using recent KLOE data concerning the gluonic content of the rj' , we are able to predict a mass formula for 
a pseudoscalar state that is pure glueball at 87%. That state would have a mass around 2.9 GeV and a decay width 
around 300 MeV if we use the pure glueball mass computed in the present pure-glue model. But large sea-quark 
effects might play an important role in that sector to decrease the mass and maybe lead to identify the ?7(1405) with 
the physical pseudoscalar glueball, as suggested in other recent works. The inclusion of sea-quark effects in effective 
models such as ours is actually a crucial issue for a better comparison to experimental results. 

We have also addressed the problem of the glueball-Pomeron identification. Our results suggest that the lowest-lying 
(even-J)++ two-gluon glueballs should no be located on the leading Pomeron trajectory, but rather on a daughter 
trajectory. The leading trajectory could rather be associated to states with a nonzero quark-antiquark component; 
such a proposal demands further and careful studies of mixing effects at the nonperturbative level to be checked. 

Finally, we have proposed a way to include finite-temperature effects in our approach, reasonably assuming that 
thermal effects cause the string tension to decrease with T. Eventually the fiux tube disappears at the critical 
temperature. Currently, lattice studies suggest that the scalar glueball mass remains constant between T = and Tc, 
while its decay width increases. The relativistic fiux tube model seems to be able to reproduce qualitatively such a 
constant mass provided that the constituent gluons gain a thermal mass which increases with the temperature. An 
increase of the thermal decay width is not excluded a priori by our approach, but further studies are needed to clarify 
the situation. Following that the thermal gluon mass is state-dependent or not, we predict the mass of the other 
glueball states to be respectively constant or significantly decreasing when approaching T^, such a prediction could 
be checked in future lattice studies. 
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masses has been plotted in gray. 
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